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The large scale interferometric gravitational wave detectors consist of Fabry-Perot cavities oper- 
ating at very high powers ranging from tens of kW to MW for next generations. The high powers 
may result in several nonlinear effects which would affect the performance of the detector. In this 
paper, we investigate the effects of radiation pressure, which tend to displace the mirrors from their 
resonant position resulting in the detuning of the cavity. We observe a remarkable effect, namely, 
that the freely hanging mirrors gain energy continuously and swing with increasing amplitude. It 
is found that the 'time delay', that is, the time taken for the field to adjust to its instantaneous 
equilibrium value, when the mirrors are in motion, is responsible for this effect. This effect is likely 
to be important in the optimal operation of the full-scale interferometers such as VIRGO and LIGO. 

04.80.Nn Gravitational wave detectors and experiments, 42.65.Sf Dynamics of nonlinear optical 
systems, 42. 60. Da Resonators, cavities, amplifiers, arrays, and rings 

I. INTRODUCTION 

The general theory of relativity predicts the existence of gravitational waves. Since gravity couples very weakly 
to matter, highly sensitive detectors are required to detect gravitational waves. Over the next decade several large- 
scale interferometric gravitational wave detectors will come on-line. These include the LIGO, composed of two 
interferometric detectors situated in the United States each with baselines of 4 km, VIRGO, an Italian/French project 
located near Pisa with a baseline of 3 km, GEO600, a British/German interferometer under construction near Hannover 
with a baseline of 600 m, TAMA in Japan, a medium-scale laser interferometer with a baseline of 300 m and with 
funding approval AIGO500, the proposed 500 m project sponsored by ACIGA |l]-||]. The large scale interferometers 
will use Fabry-Perot cavities and the ground based detectors will have arm lengths of few kilometers. There are several 
noise sources which plague the detector. Amongst them, the photon shot noise is dominant at high frequencies. It 
is reduced by increasing the amount of power of the laser source, as the noise is inversely proportional to the square 
root of the power. Therefore the cavities envisaged will operate with very high powers in their arms, tens of kilo Watts 
for initial detectors and perhaps powers as high as megaWatts in planned advanced detectors. The high power stored 
in the cavities can generate a number of nonlinear effects which would adversely affect the operation of the optical 
cavity. Here, we look into one such effect, namely, the dynamics of mirrors under the radiation pressure force. In 
earlier literature, we and others had investigated the thermo-elastic deformation of the mirrors due to the absorption 
of the power in the coatings and performed a longitudinal analysis of the cavity ||^[lO|]. Following these investigations, 
we studied the effects of radiation pressure in the cavity, in the regime when the displacement of the mirror is small 
compared with the line-width of the cavity, the cavity is servoed at resonance with a realistic servo control and the 
variation in the radiation pressure force is linearly dependent on the displacement pT[ | . The radiation pressure effects 
have also been investigated in earlier literature Jl^-p^. Here, however, since we now have a reasonably good idea 
about the instrumental parameters to be used in the large scale detectors, we expect that our analysis here will be 
important to the experimentalists. We assume that the mirrors are hanging 'freely' (there is no active servo control) 
and the radiation pressure exerts force on them which displaces them from resonance. It is very important for people 
involved in the experiments being built to have a quantitative idea of the magnitude of these effects. For instance, 
it would be necessary to take into account these effects and modify the servo transfer functions to be used for the 
control of the experiments. 

The main result of this paper is to establish that the freely hanging mirrors continuously gain energy and swing with 
ever increasing amplitude when subjected to radiation pressure force arising from the light field. This result seems 
to be in opposition with the results of Meystre et al., who found in a mirror confinement due to the radiation 
pressure. The reason for this difference of behaviour is the 'time delay' effect which is also examined in detail in this 
paper. The time delay is large for the long (kilometric) cavities, primarily studied here, while it can be neglected for 
short cavities such as the one studied by Meystre et al. and those used in most of labs. Closed form expressions 
have been given which facilitate in understanding the physics of the phenomenon. 
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The paper is organized as follows: 
In section II, we set up the equation of motion of the free mirrors. We examine the motion of the mirrors with 
the two forces (i) the radiation pressure force, (m) the force of gravity. In section III, we numerically integrate the 
equations of motion using the so called 'Phase Space Method'. We present the results for the particular case, when 
the mirrors are in the resonance positions and the laser is switched on. We observe that the amplitude of the motion 
of the mirrors increases with time and energy is pumped into the system. For very large times, when the amplitude 
of the system is also very large, so that the mirrors cross several of the Fabry-Perot resonances in one cycle of the 
pendulum, the motion approximates to that of an anti-damped harmonic oscillator. We give also, for comparison, 
the numerical results obtained for a short cavity : no anti-damping is exhibited in this case, in agreement with earlier 
results [ p^ . In section IV, we obtain analytically, under the quasi-static approximation, the phase space trajectories 
of the motion of the individual mirrors as well as the motion in the differential mode and the common mode (the 
centre of mass mode). The analytical results match with the numerical ones remarkably. In section V, we give a 
quantitative description of the phenomenon under the assumption that the velocity of the mirror does not change 
very much on the time scale of the storage time of the cavity. We find that the gain in energy is due to a differential 
radiation pressure force arising from the asymmetry, according as the mirrors are approaching each other or moving 
away from each other. We obtain approximate analytical expressions for the differential force, the time 'delay' and 
then proceed to compute the gain in energy per cycle when the mirrors encounter a single resonance or cross several 
resonances. For large amplitudes, when the mirrors cross several resonances, the system behaves like an anti-damped 
harmonic oscillator. We can then associate an effective negative Q-factor for the system. We show that the Q-factor 
depends on the input power, the finesse of the cavity and the round trip time of the cavity. Finally in section VI, we 
study the behaviour of the system when it is initially in equilibrium and goes out of lock. This may happen when the 
servo loop is suddenly opened. 



II. OPTICAL AND MECHANICAL EQUATIONS 



We consider only 'free' mirrors meaning that no servo control loop is used. The only forces acting on the mirrors 
are the radiation pressure force and gravity which manifests itself as the restoring force of the pendulum. We consider 
a single cavity with mirrors Ml and A/2 which are suspended as shown in fig.l. The input beam A enters the cavity 
from mirror All and bounces back and forth between the two mirrors. After several round trips, whose number is of 
the order of the finesse of the cavity, the field builds up inside the cavity. The magnitude of the field depends on the 
finesse of the cavity, the input power and the detuning of the cavity. The field or the power produces the radiation 
pressure force which pushes on the mirrors, driving them apart, thus changing the distance between the two. This 
in turn changes the power inside the cavity. For instance, if the mirrors were hanging in a position of resonance, the 
radiation pressure force drives the cavity out of resonance, reducing the radiation pressure force. The mirrors start 
swinging with radiation pressure force adjusting to the continuously varying length of the cavity. It is found that the 
radiation pressure force does not adjust instantaneously to the new length but lags behind the expected static force 
given by the Fabry-Perot curve (the Airy function) by a time-lag comparable to the storage time of the cavity. The 
time-lag has been called 'time delay' in earlier literature . 

The slowly varying amplitude of the field inside the cavity at time t, denoted by B{t) satisfies the following equation, 

B(t) =tiAex.p[ikxi{t)]+ RB{t-T)exp[ikL{t)], (2.1) 

where, Xi{t), Ti and ti, i = 1,2 are the positions, reflectivities and transmitivitties of the mirrors Ml and M2 
respectively; R = rir2, k = 27r/A, where A is the wavelength of the laser light, r is the round trip time and 

L(t) = 2x2{t - t/2) - xi{t) -xi{t~ t). (2.2) 

In case of the VIRGO cavity, the arm length Lq is 3 km, the round trip time r = 2Lq/c ~ 2 x 10~^ seconds, where 
c is the speed of light. This equation provides an iterative relation between the field amplitude at time t to the field 
amplitude at time t ~ t and the positions of the mirrors. We investigate the following two situations: 

1. The mirrors are hanging in the positions of resonance and the laser is switched on at the time t = tQ. 

2. The mirrors are hanging in an equilibrium state with the radiation pressure force balancing the restoring force 
of the suspension. 
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Situation 1 represents a possible experiment. When the servo-control is not operating, the mirrors will be freely in 
motion. Then, when the laser is switched on, the radiation pressure will affect the motion of the mirrors. However, 
here, we mainly deal with the case (as given in 1), when the mirrors are initially at rest and in the resonance position. 
However, many of our analytical formulae, for example, that of 'time-delay' apply to more general situations. 
Situation 2 describes the case when the interferometer is already in operation. Now if the servo-control loop is 
suddenly opened, the system will tend to become unstable. This case is also investigated. 

The equations of motion for the mirrors correspond to forced harmonic oscillator with the forcing term arising 
from the radiation pressure force. We first compute the radiation pressure forces on each mirror. For mirror Ml, the 
radiation pressure force comprises of two terms, electric field due the input laser beam, A and the intra-cavity field 
B as shown in fig.l. The radiation pressure force on Ml is, 

F,{t) = ~[R^P{t-T)-rlPo], (2.3) 

where, P{t) — \B{t)\'^, Pq = and c is the speed of light, rf is the fraction of average number of photons reflected 

by i th mirror for i = 1,2 respectively. 

For M2, the radiation pressure force is given by 

f,(,_I)=Mp„_.,. ,..4) 

The equations of motion for the mirrors with the masses, natural frequencies and damping constants, rrii, Ui, Tf, 
i = 1,2, respectively, are 

= Fi{t), (2.5) 

where x jo is the initial position of the mirror such that the separation between the mirrors before switching on the laser 

is Lq = X20 — xiQ. The full system of equations to be evolved in time are the equations from (1) to (5) (non- linearly 
coupled equations). In section III, we first carry out the task numerically and in later sections, after we have gained 
sufficient physical insight into the problem, we shall present the semi-analytical results. 
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III. THE NUMERICAL SOLUTION 



For the numerical calculation, we consider the VIRGO parameters for the suspension and the optical cavity. We 
assume mi = m2 = m ~ 28 kg., wi = W2 = i^' — 3.75 rad./sec. which corresponds to a resonant frequency of about 
0.6 Hz. The Q factor for the suspension is typically of the order of 10^. The optical parameters are r ~ 2 x 10~^ sec, 
the wavelength of the carrier wave A ~ 1.064/iim, n ~ 0.94, r2 — 1. The wave number isk = 2-k/\ and R ~ 0.94. We 
examine the behaviour of the system with the input power varying between 1 kW to 30 kW. Initial detectors will be 
operated at input powers ~ 1 kW and advanced detectors at powers of few hundred kW or even upto a MW. Note 
that these are the input powers for the main cavities after power recycling has been implemented. 

With the above values for the parameters, we find that the instabilities set in, on the time-scales of few seconds to 
few hundred seconds. Since the Q-factor of the pendulum suspension is so large, the damping in the oscillations can 
be neglected for the numerical integrations carried over the time intervals ^ % ^ 10^' seconds. Most of our numerical 
integrations range from few seconds to at most few thousands of seconds. Neglecting damping, the equations of 
motion for the mirrors become, 

Xi+uj'^{xi-Xio)='^ = fi{t), i=l,2. (3.1) 
m 

The equations of motion of the mirrors are non-linear differential equations and cannot be solved by simple methods. 

We integrate the equations by the so called 'Phase Space method' described below: 

Let A be the time-step of integration. The natural time step we assume is A = r. We assume that the forcing 
term is a constant during each time step. This assumption is not unrealistic because the round trip time interval 
for the VIRGO cavity is of the order of 10~^ seconds; note that this approximation will hold a fortiori for more 
common (much shorter) cavities. The equations simplify enormously under this assumption. Thus we can integrate 
the equations exactly within this time interval. 
The evolution of the equations goes as follows: 



3 



A . A /n(l — COSWA) ,„ „, 

a;„+i = Xn coswA + p„ smwA ^ (3.2) 

^ . fn sinwA 
p„+i = -x„sina;A + p„cosa;AH 5 , (3.3) 

where we have dropped the indices 1,2 for simpUcity and p = x/lu. Here, x represents the displaeement from the 
mean position xq. The index n represents the value of the variable at the time nA, i.e. for example, Xn = x{nA). 
The optical component has the following iterative evolution: 

mc 
2r|P„_i 

/2n = , (3.5) 

mc 

Pn = \Bnf, (3.6) 

Bn = tiAexp{ikxin) + RBn-i e^p{ikLn), (3.7) 

L„ = 2.X2„ - xi„ - xi,„_i. (3.8) 

This scheme solves the system of equations. Since it is the length of the cavity that actually matters for this problem, 
we define the variable, 

^(t) ^ k[x2{t) - xi{t)] ~ kLa, (3.9) 

and present the results in phase space plots tp/co vs ip. tp{t) is called the differential mode. 

We first consider the case when the mirrors are hanging in the resonance position and the laser is switched on. As 
the power builds up, the radiation pressure acts on the mirrors, driving them apart resulting in the detuning of the 
cavity. This reduces the radiation pressure force and the mirrors swing back. The motion is oscillatory and as we 
note, the oscillations increase in amplitude. We employ two values of input power namely, 1 kW (initial VIRGO) and 
30 kW corresponding to advanced detectors. The results are presented by the phase space trajectories of the mirrors. 
We consider four variables for the purpose, Xi, X2, tp and (j). 4>, defined later in the text in eq. (24), section IV - B, is 
called the common mode. For the range of powers considered 1 kW to 30 kW, the phase space curves obtained from 
numerical simulations, for few tens of seconds are qualitatively the same. 
We make the following general observations about the features: 

1. We observe from figs. 4 and 5 that the radius of the phase space curve increases with every cycle which indicates 
that the mirrors continuously gain energy from the input laser beam implying that the system is nonconservative. 

We also observe that the gain per cycle is not constant but some sort of a periodic function of the radius of the 
phase space diagram. We shall consider this phenomenon in detail in section V. 

2. The (static) radiation pressure force peaks when the cavity is in resonance and drops down to zero when it is 
out of resonance, (see fig. 3). The full width at half maximum (FWHM) of the radiation pressure force Flt/j) is 
about 0.06 rad, corresponding to the FWHM of the cavity resonance curve for a finesse of 50. Hence, for the 
initial stretch of the phase space trajectory, the mirrors experience the radiation pressure force whereas during 
the rest of the time, they only experience the restoring force (and the force due to the input power for Ml). The 
phase-space trajectory is circular during the restoring force regime and is deformed away from the circularity 
when the mirrors encounter the appreciable amount of radiation pressure near resonance (see fig. 4) . 

3. Since the laser is beamed in the positive x-direction, there is an asymmetry about the origin. This shifts the 
centre of mass trajectory to the positive side of the x axis (see fig. 4). The period of oscillation of the common 
mode is twice that of the period of the differential mode of the system (see fig.2). 
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4. If we let the laser beam pump in energy for large amounts of time, the amplitude also becomes large and the 
mirrors sweep across several resonances. The phase space trajectory then tends to become more and more 
circular and the motion approximates to that of a simple harmonic motion. This feature can be observed in all 
the four modes. The circularity of the trajectory implies that the motion is almost 'free'. The radiation pressure 
force has little effect because the mirrors sweep too quickly across the resonances for it to affect their motion. 
However, as we shall see that the steady gain in energy still persists. Figure 5 depicts this phenomenon. 

5. The amount of energy imparted to mirror 2 by the laser beam after getting reflected, is more than to mirror 1. 
Thus M2 swings with larger amplitude as compared to Ml (see fig. 4) (the radiation force is larger on M2). 

6. Let's turn now to the dynamics of a short cavity. The fig. 6 shows the space of phase trajectory for the differential 
mode ip, for a short cavity of length Lq — 30 cm (instead of 3 km for the other examples) and with the same 
optical and mechanical parameters and initial conditions as in fig. 5 obtained with the 3 km long cavity. The 
round-trip time is here r ~ 2 x 10~^ s, and the simulation lasts for 500 s (corresponding to a huge number of 
round trips in the cavity). At the contrary of the kilometric cavity, see fig. 5, no energy gain is found. This 
is in accord with the results of Meystre et al. , who have previously developped a theory of radiation-pressure 
driven cavities . The physical reason for this difference of behaviour between a short (common) cavity and 
a kilometric one is precisely due to the difference of length, as explained in the following sections. Time delay 
effects are thus of utmost importance for very long cavities, while they can be neglected for usual ones (and so 
have been neglected in Meystre's theory Jl2t). 



IV. PHASE SPACE TRAJECTORIES FOR THE FIRST CYCLE 



In this section, we obtain approximate closed form expressions for the equations of the phase space trajectories for 
about a period of one cycle. For the VIRGO case, this turns out to be between one or two seconds. This analysis 
could be useful in the context of the initial locking of the cavity. 



A. The differential mode ^ 



The equations of motion of the individual mirrors, equations (3) to (5) allow us to write the equation of motion of 
the system of mirrors in the differential mode as, 

4' + iu^i' = -[F2{t) - Fi{t)]. (4.1) 
m 

In the quasi-static approximation, 

p 

^ max 

l + (Wsin'V' 



Pit -r)^ Pit) (4.2) 



where = vrvT?^/ (1 — R) is the finesse of the cavity and 



_ tlPp 

^max — _ ■ \^-'^) 



Thus the equation for ■0 is. 



l + (2^)2sin^V 

where 



and 



_ 2fc(r|+^ 

i'O — 5 i^max, i^.Oj 



V'o = —■ (4-6 

mi , 1-^ 
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For the VIRGO cavity, wc have the foUowing numerical values for the quantities: 
V'o - 0.88Po/10kW, J'-bO and Fq ~ 62.6Po/U)kW. 

For ^ ^ 1, the term on the right hand side of equation (18) is non-zero only when tfj 1. With the approximation 
sin^ ip, we can easily integrate equation (18) to get the phase space trajectory of the mirrors, 

^ + (V; + ^of = ^ tan-i(^^) + (4.7) 

For low powers like 1 kW, the approximation sini/; ~ works remarkably well and agrees with the numerically 
obtained phase space trajectory. In order to compare the analytical and numerical results, we compare the maximum 
value of tjj, namely V'maxi of the trajectories for various input powers in fig.7. 

1. When ^ is small that is near resonance, the equation of the trajectory reduces to 

^ = 2i?bV- (4.8) 

Thus the trajectory is parabolic in shape and passes through the origin. 

2. When ip ^ 1, the trajectory is a circle as expected since there is hardly any radiation pressure force acting on 
the mirrors. The equation of the trajectory in this regime is 



^ + (V + Vo)' = ^+Vo- (4.9) 



When the input power is very large 50 kW, the trajectory does not maintain this simple shape. For example, 
when Po = 50 kW, the trajectory is as shown in fig. 8. 

This is because the high powers make the mirrors cross several resonances in the first cycle itself and consequently 
the trajectory has more complex behaviour. We do not pursue this case here. 



B. The common mode (j> 

In this sub-section, we study the motion of the center of mass of the system of two mirrors. We define the center 
of mass coordinate of the two mirrors as, 

(p = k{xi - Xio) + k{x2 - X2q)- (4-10) 

The equation of motion of the system in the center of mass coordinate is, 

:j> + u;^^=^[Fiit)+F2{t)], (4.11) 
where for Fi{t) and F2{t) are given by equations (3) and (4). The equation of motion takes the form. 



1 + (^)^sm V 



,2 



where 



2fc(ri - R')tlPQ 



(4.13) 



For VIRGO parameters, Fc — ^-^ ^okw "^^^ above equation is a second order differential equation and is coupled 
nontrivially to the tp mode. The strategy we adopt is to study the motion of the center of mass in different regimes; 
(1) near the resonance and (2) away from the resonance. The full trajectory is obtained by matching the solution in 
the region of the overlap. 
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1. For < < 0.5; sin^ ip, equation (26) becomes 



Fc 

1 + 



V'o) = ^^^. (4.14) 



Further for -0 < 3 x 10 ^, {^^)^ ^ 1; we neglect (2^)^ as compared to 1 and obtain, 

^+</.2 = 2(Fe + Vo)</'. (4.15) 

For the VIRGO cavity and input power 1 kW, Fc ~ 0.4. Neglecting the quadratic term in (j) as compared to the 
linear term, we see that the motion of the center of mass describes a parabola for low values of tjj and (j), 

^ ~ 0.970. (4.16) 

We compare the slopes of the phase space diagrams in the differential mode as well as the common mode. The 
phase space curve is more steeper in the common mode as compared to the differential mode. 

2. Away from resonance, tp > 0.5, the power stored inside the cavity is almost zero. The equation of motion of the 
center of mass is given by, 



= V'o- (4.17) 



The solution for the initial conditions (j) = 4>o,4> = 4>o is, 



4 + (<^-V'o)' = 4 + (<^o-Vo)'. (4.18) 

This solution must be matched to the solution in case 1. Suppose we match the solution at 4>q = 0.5, from equation 
(29) we have ^/w ~ 0.7. This gives the approximate solution. In general, the solution in this region is, 

^ + {.4>- V'o)' = 2(Fe + Vo)<^o + (<^o - ^of - 4>l- (4.19) 
The equations (29) and (33) describe the full solution for this mode. 

The phase space trajectory for the individual mirrors can be obtained from the motion of the mirrors in the 
differential and the common modes. The trajectories of the xi and modes are shown in figures 2 and 4. 



V. ENERGY CONSIDERATIONS AND ANTI-DAMPING 



In the previous section, we examined the motion of the mirrors under the radiation pressure force. We noted that 
each time when the mirrors cross a resonance they experience an impulsive force. In this section, we analyze in detail 
as to how the system gains energy as each resonance is encountered. 



A. Quasi-static approximation 

We start with the quasi-static approximation in which the system is conservative, thus there is no net gain in energy. 
We then phenomenologically introduce the 'time delay', Tiag which now leads to gain in energy. We thus obtain Tiag 
in terms of the cavity parameters and this gives us an equation for an anti-damped harmonic oscillator. The energy 
gain can be obtained as shown in the following sections. Finally, for large times, the gain can be expressed through a 
negative Q of a harmonic oscillator. In the quasi-static approximation, we assume that the mirrors arc moving 'slowly' 
that is the intra-cavity power has time to adjust itself to the slowly changing positions of mirrors, i.e. Tip — R. 
Thus neglecting r from the equation of motion we obtain, 
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^ + i^ + i;o)^Fs{t^), (5.1) 

where tpo as given in the previous section gives the constant displacement due to the constant input power Pq from 
the laser beam and Fs{tp) is the radiation pressure force in the static case is given by, 



l + (^)2sin^^' 

where Fq is given in equation (19). The dimensionless energy of the system is an integral of motion and is obtained 



E=-^ + VsusW + Voii^) + VradW, (5.3) 



where the potentials are given by, 

1 



Vsus{^) = (5.4) 



Vo[^) = V'oV', (5.5) 

K-ad(v-) =- r Fsi^P)d^. (5.6) 
Jo 

The system is conservative in this approximation. These results have been obtained in the earlier literature ]l^ . 

However, we note from the results of numerical simulations that the system indeed gains energy with every cycle of 
oscillation. Moreover, this gain occurs when the radiation pressure force is appreciable i.e. when the system is near 
resonance. Most importantly, we observe from the numerical simulations that radiation pressure force lags behind its 
quasi-static value by a 'time-lag' which we denote by Tiag- We find that in the case of VIRGO, Tiag varies from about 
16t to 30r as one climbs up the resonance curve from half its maximum to the maximum. Tiag is of the order of the 
storage time of the cavity. With this in mind we write the force F{t) ^ Fg{t ~ Tiag) and obtain the following equation 
of motion, 

^ + {^ + i^0) = Fsm-Tlag)). (5.7) 

Taylor expanding the forcing term to the first order we obtain, 



, _ dF, 



2 +^'-9-rfV' + (V' + ^o) = i^.(V'), (5.8) 



as the equation of motion for the system. The ip term in equation (41) is responsible to the gain/loss of energy of the 
system. This solely depends on the sign of since Tiag is always positive. If we start the system from resonance 

at ip — 0, ^ starts increasing slowly and in the region, < and the system experiences more force than what 
it would have been in the quasi-static case. This excess force is manifested in an excess amount of energy AE. In 
general, the energy gained/lost over a certain amount of time is given by, 

AE=- riagW-T^^'d^^ = - Tiag-r^^P^dt. (5.9) 

It is clear from equation (42) that when < 0, the energy gain AE > and vice-versa. But since the radiation 
pressure force always gives a kick in the positive -ip direction, the is always larger when < 0. In other words, 
the amount of energy gained by the system is more than the amount of energy lost. Hence, there is always a net gain 
in energy when the system crosses a resonance. We need to compute this extra force AF which gives rise to the gain 
in the energy to the first order in Tiag- We obtain. 



AF=-^^Tiag. (5.10) 
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We note that AF depends on the phase velocity and Since ~ when the system is away from resonance, 
AF comes into play only in the region of resonance. Another way of describing the fundamental asymmetry is to say 
that for the moving mirror, the wavelength of the laser light is modified by the Doppler effect. When the mirrors are 
approaching each other, the apparent frequency of the light for the cavity is increased, or equivalently, the line-width 
of the cavity seen from the laboratory frame is narrower and vice versa when the mirrors move away from one another, 
the line-width is seen to be broader. The consequence is that the braking force (when the mirrors move against the 
light) acts for a shorter time than the accelerating force, when the mirrors are moving away from each other. Over 
one cycle the energy difference is positive and there is a continuous increase of mechanical energy which comes from 
the laser. 



B. The computation of the time-delay 

In this section, we obtain a closed form expression for the 'time-delay', riag under the approximation that the relative 
velocity of the mirrors docs not change much during the storage time of the cavity. This is observed in the numerical 
simulations and hence the approximation may be justified. Tiag is the most important quantity for computing the 
rate of gain in energy. The intra-cavity radiation field at n-th time instant S„ is related to the intra-cavity field at 
n + 1-th time instant via an iterative relation, 

Bn+i = tiA + Rexp{2itl>n)Bn, (5.11) 

where, = B^iit), (,';„ = ^i'(nr). Wc compute the field i?„ from initial time to, which we take to be zero. We take 
the time step for the iteration to be the round trip time r. In order to compute the intra-cavity field, it is necessary to 
know the temporal behaviour of tp. For the static case, when the mirrors are stationary, we get the equilibrium field 
Bs- When the mirrors are moving, the approximation now comes into play, namely, we assume that up is constant 
over the storage time of the cavity. 

The Taylor expansion of V'fe to the first order around t = Q is, 

-ijjk = ipo + kTijjo, (5.12) 
where ipo = V'(^ = 0)- Iterating equation (44) n times starting from t = when B = Bq, we obtain, 

Bn = hA i?"e^'^^=— '^'^ + ii^e^'Sr^o ^'=Bo. (5.13) 

m=0 

For large n, the last term in equation (46) tends to zero and the electric field amplitude at the n-th iteration is given 
by, 

Bn = tiA BJ^e^^2Zl=n-^'l'\ (5.14) 

m=0 

• If the relative separation between the mirrors is constant in time, ipk = i>o, we retrieve the equilibrium field Bgi 

OO J. A 

B, = t,A Y R-^e^'-'^'o = (5.15) 

m=0 

• Turning now to the case in hand, when iph varies linearly with time, the sum in the exponential for large n is 
approximately given by 

n—l ^ 

^ ■^fc ~ mipn - -^m^TipQ. (5.16) 

k=n—m 

We combine equations (47) and (49) to obtain, 

B„ = tiA Y Ji'"e2™'^"e-™'^'''«. (5.17) 

TO=0 
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If ip is small then the exponent can be linearized. Then it is possible to express the field as the sum of the static 
term and the remaining part AB which corresponds to the time lag. 
We write 

Bn^Bs + AB, (5.18) 

where 

n-l 

AB = -iTtiA4>o J2 -R"e^""'^"TO^- (5.19) 

?TI = 

Summing the arithmctico-gcometric series jl7| yields, 

AS = -2^Mr^^^^__. (5.20) 

The corresponding power AP is given by 

AP = 2i?e(i3„ABJ ^ _ ^ ■ (5-21) 

Refering to equation (19), we get the expression for the dimcnsionless extra force, 

A^.^M±^AP. (5.22) 

The above expression of power is for values of V' near resonance when -0^1. When ij; ^ mr, the same expression 
can be replaced by the ip — nir under similar approximations. The effective 'time delay' Tiag is now obtained 
from (43), (54) and (55). 



- 2ri? 

"'"^ ^ (l-i?)(l + (2^)2^2)- ^^-'^^ 



We note from the equation (56) that the effective time lag is maximum at the resonance and starts decreasing as 
one goes away from the resonance value. For VIRGO parameters, we compute the value of the Tiag. At FWHM, i.e. 
ip ~ 0.03, Tiag ~ 16t, going up to Tiag ^ 30t as one approaches resonance. 



C. Energy gain near the resonance at ^ — 

In this section, we examine the motion of the mirrors for small amplitudes. Initially before the laser is switched on, 
the cavity is at resonance i.e. the two mirrors are separated by an integral multiple of tt in phase (^^ is an integer). 
The system starts from tjj = and the motion is allowed to evolve with time. We restrict the amplitude to \ip\ < tt 
and study the system in this regime. Our goal is to compute the net gain in energy, AEcyde during one cycle of 
oscillation. We have, 

AE,y,i, = 2 AF{^)di^, (5.24) 

where 







A^(^) = ^^^^^^°tt ■ (5.25) 

We note that < in this regime. We set an arbitrary cut-off tpi (when AF <C 1) as the system moves away from 
the resonance. The factor of two is because the force is encountered twice during the cycle. 
The ip is obtained from the energy balance equation. 
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= r F,{i,)di, ~ ^Fotan-i(^). (5.26) 

Thus from (57) and (58) we get, 



AEeyde = 8wT-y/ ^i^o^'/(<^i), (5.27) 



where, 



, , /-^i (/.(tan- 1(^)1/2 ^ ^ 

and = 2Til)\/'K. The cut-off should be away from the resonance and we find for VIRGO, ~ 5 is an acceptable 
value. 

The numerical and the analytical results are compared in fig. 9 by plotting the gain in energy per cycle for various 
input powers. 

D. Energy gain for large amplitudes 

As the mirrors gain energy, they swing with ever increasing amplitude and sweep over several resonance peaks of the 
Fabry-Pcrot cavity. The system gains energy at every resonance peak and thus the energy gained per cycle is the sum 
of the energy gained at each resonance encountered. The peaks are encountered at = nir, with rimax ^ n > rimm, 
nmin < and rimax > 0. The total number of resonances encountered by the mirrors is Umax + \nmin \ + 1- Also it 
is observed that Umax > I'^mml due to the built-in asymmetry arising due the laser power pumped in the positive 
a;-direction (see fig. 3). 

Let be the energy gain at the n-th resonance, then the total energy gained per cycle AEcyde is given by, 

^E^ycie = ^ Ai;„, (5.29) 

where 

A£;„ = AEn+ + AS„- , (5.30) 
and i^En-, A£„+ is the energy lost or gained respectively and given by the following expressions, 

/^En+ = I AFdTp > 0, (5.31) 

J nir 
rnTT 

AE„_ = / AFdtp < 0. (5.32) 



Combining equations (57,58,62 to 65) we get the energy gain at nn as, 

" Tr^(^"+ - i'^ (i + (L).)3 ^^' (5.33) 

where ipn- and 'ipn+ arc the relative phase velocities of the mirrors for tp < mr and tp > mr respectively. The energy 
gained at different resonance positions of the mirrors is different because the ip is different at different resonances, tp 
is maximum when |n| is small and becomes small when n approaches nmin or Umax- Thus the net energy gain per 
cycle is 

l^E^ycle ^ )' E (^"+ - (5-34) 
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The next task is to compute the {ipn+ — ipn-), the increase in phase space velocity while crossing the resonance 
position as a function of n. For large values of n, the radiation pressure effect reduces remarkably as is observed in 
fig. 3. Hence wc take the static force equation to compute the increase in the phase space velocities. We integrate 
equation (41) neglecting the anti-damping term obtaining the change in kinetic energy as, 

1 • ■ f°° I 

^{^1+ ^l.) = Fo - 2V^iV'o - 2n7rVi. (5.35) 

Equation (68) reduces to 

(^„+ - V>„-) ^ ^-4V;i(V>o + n7r) 
w 2i){n'K)/uj 

where we have approximated V'n+ + ^/'n- — 2tp{mT). 

To compute '^{n'K)/oj analytically, we consider the dimensionless instantaneous energy as given by equation (36), 

where J(V') = {i-ry'+arsi-d:^ ^ ' ^® approximate the integral for large motions as J{tp) ^ ^"^2 ~ for 
n crossings of the resonances. Assuming r2 = 1 and ri = R, we have, 

f^ + {^-i,,f = {2E + i,l), (5.38) 

where Vc = mu^c ^ wkw (^^'^ VIRGO parameters). The phase space trajectory is a circle centred around Vc with 
radius 

p = v/2¥TW- (5-39) 

Equation (71) gives tisul as 



We rewrite the gain in energy per cycle as, 



Vp2-(n7r-Vc)2. (5.40) 



AEcycle ^ j^)3 \'^) ^\^max,nmin), (5-41) 

where Umax is the greatest integer not greater then (f/^c + p)/iT and rimin is the smallest integer not smaller then 
(V'c — p)/7r. We have, 

Umax 

S{nmax,nmin) = ^ ginir), (5.42) 

Tl=nm.in 

where 

^-2V>i(^o+n^) 
- (riTT - ipcY 

For large motion of the mirrors, since the mirrors cross many resonances, we may replace the sum by an integral over 
n. Changing over to the variable M = n-K — ipc, in which the system appears more symmetric about the origin, we 
have, 



dM (5.44) 
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where a = - 2±L(^g + and /3 = Mmax and M„j„ correspond to Umax and rimin respectively. 
We note that p satisfies the following inequalities, 

M^ax< P<M^ax + l. |M™„| <p< |M™„| + 1. (5.45) 

We observe that when n is large, the difference between |Mmi„| and Mmax is small, of the order of 1 or 2 times tt. 
We denote the difference, by 5M where, 

5M = Mmax - \Mmin\ = Mmax + Mmin- (5.46) 

The integral in equation (77) splits into three parts 5 = S'l — 5*2 — 6*3, where 

Si =2asm-^ Mmax/ P, (5.47) 

2a 

Vr ^^-^niax) 

The dominant term is which gives the general behaviour and shape of the curve shown in fig. 10. While ^2 produces 
a small kink in the curve, the effect of ^3 can essentially be ignored. In case of the VIRGO cavity and Pq = 30 kW 
and p ~ 20, Si ~ 16.7, ^2 =± 1.75 and ^3 =± 0.26. For the initial VIRGO detector, Pq = 1 kW, and p ~ 4.0, 5i ~ 0.56 
while S2,Sic:t 0. 

Considering only Si, the energy gain per cycle is approximately given by, 

A 7-^ ^^^max . —1 Mmax n ^ ^ ^ a r , /r rr\\ 

AEcycle = Sm , Mmax < P< Mmax + 1, (5.50) 

TT p 

where, 

AEmax = _ j^)3^2 " " ('^O + V'c)J • (5-51) 

We observe the following features in the profile of AE cycle' 

• At the resonance position , p = Mmax the energy gain is maximum and equal to AEmax- For the VIRGO cavity 
specifications and input power of 30 kW, AEmax — 3.8. When Pq = 1 kW, AEmax — 4.2 x 10~^. 

• Equation (83) shows that the energy gain per cycle is a decreasing function of p. The energy gain decreases till 
the next resonance is crossed, where it suddenly increases to AEmax- As the amplitude increases, p increases 
from Mmax to Mmax + 1) the mirrors sweep across the resonances a little faster which deprives them from 
gaining the full energy AEmax- ^Ecyde therefore reduces from AEmax to AEmin, where 

A ri 2AEmax . —1 Mmax /r rr»\ 

AEmin = sm'— — -. (5.52) 

TT Mmax + J- 

• The energy profile in this range of p i.e. Mmax to Mmax + 1 can be approximately given by. 



{AEmax - AEcyclef = ^^^P^iP ' Mmax)- (5.53) 



We note that the minimum value is, 

AEmin ~ AEmax ( 1 - ( ) ) , (5.54) 

V \-K^Mmax'' ^ 

and it tends to the maximum value AEmax as Mmax becomes very large. 

• In fig. 10, it is seen that there is a kink after AEmin is reached . The kink occurs because Mmin reduces by 1 
when the mirrors cross yet another resonance on the negative side. This is accounted for by the second term. 

Since the E p^, the rate of increase of p per cycle, denoted by Ap is given by Ap = AR, 
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E. The negative Q-factor 



In the previous section, we have seen that for large values of the amount of energy gained by the system of two 
mirrors per cycle is a constant and tends to AEmax- The system is an anti-dampcd harmonic oscillator which gains 
on an average, constant amount of energy per cycle. We may therefore associate a negative quality factor —Q, where 
Q > 0, with the system which describes the anti-damping. In this section, we endeavor to study the behaviour of Q. 
The equation of motion of an anti-damped harmonic oscillator is given by, 

^-^+X = 0, (5.55) 

where Q is constant. The solution for x is of the form 

X ~ Xoe'^*/^'5e±''"*. (5.56) 

where x — ''P ~ i^c- The amplitude is p ^ Xoe"*^^*^ when p is sufficiently large so that the phase space trajectory is 
approximately circular and can be compared to a simple harmonic oscillator. Since a constant amount of energy is 
gained per cycle, the Q will be a function of time. However, we can still describe the system by an average Q taken 
over a cycle, which we denote hy < Q >cycie, since the change in Q during one cycle is small. Moreover, we assume 
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'''<Q>^y':i' ^ \ i c the fractional variation of < Q >cycie over a period of a cycle can be ignored. Thus we 



obtain a W.K.B. solution for p{t) as, 

p(t)=p(to)exp - / — — , (5.57) 

'-■^ J to ^ ^cycle-' 

where to is some fixed but arbitrary initial time instant and /o(io) is the radius of the circular phase-space trajectory 
at to- We have also assumed that < Q >cycie^ 1- From equations (72) and (90) the energy of the harmonic oscillator 
is given by, 

E{t) - E{to) = lp\to) [eM[ ^ Q^^^J - l] • (5-58) 
Since the energy gain is a constant and during a cycle equal to AEmax, we can equate to AEmax to obtain, 

E{t) - E{to) = "^^^^{t - to). (5.59) 
The time evolution of p is obtained, which yields, 

Pit) = p{to) [1 + ^%^a;(t - to)] (5.60) 
Equating the logarithmic derivatives of (90) and (93) we obtain, 

< Q >cvcle (t) =< Q >cycle (to) + w(i - to), (5.61) 

where < Q >cycie (to) = " ■ We observe that both the energy and the < Q >cycie increase linearly with time 

while the amplitude p increases as t^^'^. < Q > cycle (t) depends through AE^ax on the input power, finesse and the 
round trip time. 

For Pq = 30 kW of input power, the trajectory more or less obtains a circular shape when p{to) ^ 15. For the 
VIRGO parameters, 

AEmax ~ 0-42(y^)' ~ 3.79. (5.62) 

Thus < Q >cycie (to) ~ 186 and so 

< Q >cycle (t) ~ 186 + L0{t - to). (5.63) 
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Further, if we also consider the effect of the damping of the suspension then the hmit cycle will be approached 
when < Q >cycie'^ Qsus ~ 10^ for VIRGO. < Q >cycie will attain this value after uit ~ 10® which corresponds to 
little more than 3 days. The corresponding amplitude is given by, 

< Q > cycle [to) 

This simple analysis will have to be modified when the limit cycle is almost reached, that is when, < Q >cycie^ Qsus- 
Here, however our goal was to estimate the time it takes to reach this stage. The above analysis is then adequate for 
the purpose. However, when the servo operates or otherwise such a situation is unlikely to arise because other effects 
such as mirror tilting etc. will be important long before and then the dynamics will be completely different. 



VI. THE EQUILIBRIUM CASE 

Lastly, we examine the case when the mirrors are in equilibrium under the radiation pressure force and the restoring 
force. We have to determine whether the equilibrium is stable or unstable. To this end, we perturb the equation (41) 
about a given equilibrium point. Writing ijj = ip^q + Sip and linearizing, we get, 

5^) - —Sip + fl^6tp = 0, (6.1) 



' eq 



where, 



and 



If ipeq < 0, fl^q < 0, the sigu of Teg docs not matter and the instability grows exponentially. If on the other hand 



tpeq > 0, then although > 0, Teq > and this leads to gradually growing oscillations until the pendulum tips over 
the maximum. We plot the phase space trajectory for the input power of 1 kW in fig.ll. We conclude from this that 

the cavity is always unstable when radiation pressure forces act. The time-delay plays a crucial role in making the 
system unstable. The negative Q-factor for the motion of the mirrors near the equilibrium position, is given by. 



g = --f2egTeg. (6.4) 



For an input power of 30 kW, with VIRGO cavity parameters, the equilibrium position of the mirrors near the 

resonance at zero is ipeq — 0.25. The corresponding values of the other quantities are ^l^q — 4.75 and r^q ~ 10'' 
seconds. The negative Q-factor, Q ~ —2.3 x lO**. Whereas for initial VIRGO, Pq — 1 kW, tp^q — 0-16 thus il^g ~ 1.67, 
Teg — 3.1 X 10^ seconds, the negative Q-factor, Q ~ —2.6 x 10^. 



VII. CONCLUSION 



We have analysed the effect of radiation pressure on the freely hanging mirrors (no servo loop) suspended in the 
laser interfcromctric optical cavities. After numerically evolving the full set of equations of motion with respect to 
time, we find that the amplitude of the mirror oscillations continuously increases as time progresses. We introduce 
the 'time delay', that is the time taken for field to adjust to the motion of the mirrors, in a phenomenological way 
to explain the observed gain. We conclude that the gain in energy is due to the differential radiation pressure force 
arising from the asymmetry depending upon the motion of the mirrors. From another viewpoint, we can also explain 
the gain in energy qualitatively by the Doppler effect. With respect to the mirror, the frequency of the incoming laser 
beam is higher as compared to that of the outgoing laser beam due to the Doppler effect. The deficit of the energy of 
the laser beam after getting reflected from the mirror can be looked upon as the energy gained by the mirrors. The 
values of the energy gain per cycle are computed analytically under the reasonable assumption that the mirrors are 
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not accelerated within the time scale of the storage time of the cavity. For VIRGO parameters, the analytical values 
agree remarkably with the numerical values. The interesting point to note is that the motion of the mirrors approaches 
that of an anti-damped harmonic oscillator with a constant gain in energy as time progresses which implies that the 
mirrors move too quickly to get affected by the radiation pressure force. The negative Q-factor of the anti-damped 
oscillator depends on the input power, the finesse and the round trip time of the cavity and increases linearly as 
a function of time. The analysis will have to be modified when the negative Q-factor becomes of the order of the 
damping Q-factor of the suspension fibre, if such a case can arise. 

In this paper, we have shown that the radiation pressure force makes the freely hanging mirror unstable for all 
values of the input power and irrespective of the initial conditions. 

The above analysis is relevant in the event, when the interferometer is in operation and if the servo loop is suddenly 
opened. Then the motion of the hanging mirrors can be deduced from the above analysis. This analysis will be helpful 
in designing a servo-control which can prevent this instability. In a previous work [ pl| , the servo-control was included 
in the linear regime of the Fabry-Perot curve assuming the transfer function for the servo given by Caron et al. . 
Their work sets the stage for analysing the system in the non-linear regime as well, but it is then needed to know how 
servo-control operates in the full regime. 
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FIG. 1. Schematic diagram of the cavity and the intra-cavity fields. 
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FIG. 2. Motion of the mirrors for the modes kxi, kx2, tp and ^ as a function of time for input power of 1 kW. 
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FIG. 3. Radiation pressure force profiles for the modes kxi, kx2, and (p in dimensionless units for input power of 1 kW 
for one cycle. The conversion factor for converting the dimensionless force into Newtons is mw^/fc. 
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FIG. 4. Phase space diagrams for kxi, kx2, ■0 and 4> for 1 kW of input power and integration time of 3 seconds. 
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FIG. 5. Phase space diagram for tp for the input power of 30 kW and integration time of 500 seconds. 
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FIG. 6. Phase space diagram for ^ for the input power of 30 kW, integration time of 500 seconds and for a short cavity 
(I/O ~ 30 cm). Unhke the long cavity, studied here, there is no energy gain : the phase-space trajectory remains remarkably 
stable. 
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FIG. 7. Comparison of the values of tpmax obtained analytically (smooth curve) and numerically (open circles) for input 
powers of 1 kW, 5 kW, 10 kW and 30 kW. 
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FIG. 8. Phase space diagram for the tp mode for the input power of 50 kW and integration time of 3 seconds. 
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FIG. 9. The gain in energy per cycle obtained analytically (smooth curve) and numerically (open circles) for input powers of 
1 kW, 5 kW, 10 kW and 30 kW. To convert the dimensionless energy gain in Joules we multiply by the factor of mw^/fc^ ~ 11.2 
pico- Joules. 
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FIG. 10. The dimensionless energy gain per cycle AEcycU as a function of the amplitude of the ip mode for the input power 
of 30 kW. The plot is analytical. 
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FIG. 11. The phase space diagram for the equihbrium case for the tp mode. The input power is 1 kW and the integration 
time is 4 seconds. 
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